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ABSTRACT 

A star acquires much of its mass by accreting material from a disc. Accretion is probably not 
continuous but episodic. We have developed a method to include the effects of episodic accretion in 
simulations of star formation. Episodic accretion results in bursts of radiative feedback, during which 
a protostar is very luminous, and its surrounding disc is heated and stabilised. These bursts typically 
last only a few hundred years. In contrast, the lulls between bursts may last a few thousand years; 
during these lulls the luminosity of the protostar is very low, and its disc cools and fragments. Thus, 
episodic accretion enables the formation of low-mass stars, brown dwarfs and planetary-mass objects 
by disc fragmentation. If episodic accretion is a common phenomenon among young protostars, then 
the frequency and duration of accretion bursts may be critical in determining the low-mass end of the 
stellar initial mass function. 

Subject headings: Stars: protostars, formation, low-mass, brown dwarfs - Accretion, accretion disks 
-Methods: Numerical - Hydrodynamics - Radiative Transfer 



1. INTRODUCTION 

Stars form from the collapse of gravitationally un- 
stable cores in molecular clouds. Interstellar turbu- 
lence ensures that such cores have a net angular mo- 
mentum, and so protostars form attended by accretion 
discs; much of the collapsing material first infalls onto the 
disc and then spirals into the protostar (jTerebev et al.1 
1984). If the disc cannot efficiently transport angular 
momentum outwards, accretion onto the central proto- 
star is slow, and the disc grows in mass. The disc may 
then become gravitationally unstable and fragment (jBatel 
l2009al lAttwood et al.ll2009l ). spawning secondary proto- 
stars, primarily of low-mass (iStamatellos fc Whitworth! 
l2009allStamateilos et al.ll2010f ). i.e. low-mass hvdrogen- 
burning stars (0.08 M Q to ~ 0.2 M Q ), brown 
dwarfs (0.012 M Q to 0.08 M Q ), and even planetary- 
mass objects (< 0.012 M^). This mechanism 
ijWhitworth fc Stamatellos! I2006D is the only proposed 
mechanism that addresses in detail all the critical obser- 
vational constraints, i.e. the shape of the low-mass end of 
the initial mass function, the brown dwarf desert, the bi- 
nary statistics of low-mass objects, an d the formation of 
free-floating planetary-mass objects (|Stamatellos et al.1 
I2007at IStamatellos fe Whitworth! [2"009aD . 

However, as mass flows towards the protostar, gravi- 
tational energy is transformed into thermal energy due 
to viscous dissipation in the disc and at the accre- 
tion shock around th e protostar. It has been pointed 
out (IKrumholzl [2006t iBatel I2009bt lOffner et all |200~ 



brown dwarfs form in simulations of fragmenting clouds, 
whereas with radiative feedback a smaller numb er of 
brown dwarfs form. However, the simulations of IBatel 
( 2009b) include radiative feedback due to the conversion 
of gravitational energy to thermal energy only down to 
the sink radius, i.e. 0.5 AU; thus, IBatel (|2009b[ ) simula- 
tions only account for a small fraction (~ 2%) of the total 
radiative feedback from young protostars and represent a 
lower limit on the effect of radiative feedbac k in disc frag- 
ment ation and low mass-star formation (cf. lOffner et al.1 
2009'). When the radiative feedback from young proto- 
stars is fully included then the formation of low-mass and 
brown d warfs by disc fragm entation is almost fully sup- 
pressed ([Offner et al.l 120091 their Fig. 6). Furthermore, 
based on the effects of radiative feedback, lOffner et al.1 
(2010) conclude that disc fragmentation is not a major 
contributor to low-mass stars and brown dwarfs. More- 
over, disc fragmentation is currently the only model that 
predicts the statistics of low- mass s tars and brown dwarfs 
()Stamatellos fc Whitworth! I2009al) . Hence, the forma- 
tion of low- ma ss stars (Z; 0.2 Mg), wh ich, if one includes 
brown dwar fs (iWhitworth et al.ll2007h . constitute ~60% 
of all stars (|Kroupall200li r remains an open problem. 

Previous studies have assumed that the accretion 
of material onto protost ars, and hence their radiative 

' " Batj 12009b!: 
However, 
theoretical evidence 



feedback, is continuous dKrumholz 



Urban et al.l l2010t iKrumholz et al.1 120101: lOffner et al 
20U¥ that this accretion-related energy heats and sta- 



lOffner et all 120091: IKrumholz eta! 
there is strong observational and 
that accretion of material is often episodic dHerbie 



Dopita 1978; Reipurth" 1989 : lHartmann fc Kenyan 



1977 



1996 



Greene et al. 2008; Peneva et al.ll2010t ): it is concentrated 



bilises the disc, thereby suppressing frag mentat i on and 
low- mass star formation. More specifically, IBatel J2009b) 
finds that when one ignores radiative feedback too many 
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in bursts that last from a few tens to a few hundred years. 
In particular, FU Ori-type stars exhibit sudden large (~5 
mag) increases in their brig htness, and then fade away 
over decades or c e nturies (lHartmann fc Kenvonl 119961 : 
iGreene et afl 120081 : iPeneva et all 120101 ) . The estimated 
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accretion rates during these bursts are up to ~ 5 x 
10~ 4 M© yr -1 , so each burst can deliver a few 10 -2 M© 
onto the protostar. Further evidence for episodic accre- 
tion comes from the periodical ly spaced knots seen in 
bipolar jets (e.g. IReipurth| [l989) . Bipolar jets are driven 
off along the protostellar rotation axis by the energy re- 
leased as material spirals in and accretes onto the proto- 
star; periodically spaced knots therefore imply episodic 
accretion. Finally, the luminosity problem provides indi- 
rect observational support for episodic accretion. By the 
end of the Class phase (~ 10 5 yr), half the protostar's 
final mass has been accumulated. Thus for a final mass of 
1 Mq, the mean accretion rate onto the protostar must be 
~5x 10 -6 M Q yr -1 , and the mean accretion luminosity 
must be ~ 25 L©. This is much larger than the observed 
bolometric luminosities of typical sola r-type protostars 
(|Kenvon et al.lll99tt lEvans et al.H2009H . Episodic accre- 
tion mitigates this problem, as the luminosity is only 
large intermittently, and most protostars are observed 
between bursts. 

The cause of episo dic acc retion is uncer- 
tain (IBonnell fc Bastien] Il99l IBell fc Linl Il994 ; 
Vorobyov fc Basul I2005T IZhu et all 120071 I2009allb . 
2010bllari. but t he most compelling explanation 



(|Zhu et a"H l2009aft involves the interplay between 
two different modes of angular momentum transport. 
The outer disc is sufficiently cool that it is susceptible 
to gravitational instability (GI); as well as creating 
low-mass companions in the disc, GI also results in 
gravitational torques that transport angular momentum 
outwards, thereby allowing material to spiral inwards 
towards the central protostar. In contrast, the inner 
disc is too hot for GI, but, if it becomes hot enough, 
the gas is thermally ionised to a sufficient degree that it 
couples to the magnetic field and the magneto-rotational 
instability (MM) kicks in; the MRI then transports 
angular momentum outwards, allowing matter to spiral 
in further and onto the central protostar. The upshot is 
that matter accumulates in the inner disc, and steadily 
heats up until it is sufficiently ionised to activate MRI, 
whereupon the accumulated matter spirals into the 
central protostar, giving rise to an outburst. Once most 
of the accumulated matter in the inner disc has been 
deposited onto the central protostar, the temperature of 
the remaining matter falls, the MRI shuts off, and the 
accumulation of matter resumes. 

In this paper, we present a sem i-analytic prescrip tion, 
which is based on the model of IZhu et al.l ()2010bl ). to 
incorporate episodic accretion in radiation hydrodynamic 
simulations of star formation. We perform simulations of 
the collapse of a star forming molecular cloud and show 
that episodic accretion plays a crucial role in promoting 
disc fragmentation and low-mass star formation. 

2. COMPUTATION METHOD 

We use the SPH code seren 
(Hu bber. D. A.. Battv. C. McLeod. A.. Whitworth. A. Pi 
2011]), in which the chemical and radiative processes 
that regulate the gas temperature ar e treated with the 
flux-lim ited d iffusion appro x imatio n of lSta matcllo s et all 
(|2007bft and iForgan et all (|2009D . The net radiative 
heating rate for the particle i is 



The positive term on the right hand side represents 
heating by the background radiation field, and ensures 
that the gas and dust cannot cool radiatively below 
the background radiation temperature T BGR . a SB is 
the Stefan-Boltzmann constant, is the mass-weighted 
mean column-density, and R R (pi,Ti) and K r (pi,Ti) are 
suitably adjusted Rosseland- and Planck-mean opaci- 
ties. The method takes into account compressional 
heating, viscous heating, heating by the background 
radiation field, and radiative cooli ng/heating. The 
method has been extensively tes ted (Stamatello s et all 
I2007bt iStamatellos fc Whitworth] 1200811 In particular 
it rep roduces the results of both detailed 3D simula- 
tions (iMasunaea fc Inut suka 200"ol: iBoss fc Bodenheimerl 
119791 : iWhitehouse fc Bate! 120061 ) . and als o -unlike other 
schemes - a nalytic test calculations (Spiegel 1957; 
IHubenvl ll990! The gas is assumed to be a mixture of 
hydrogen (70%) and helium (30%) , with an equation 
of state (Bla ck fc Bodenheimerl I1975I ) that accounts (i) 
for the rotational and vibrational degrees of freedom of 
molecular hydrogen, and (ii) for the different chemical 
states of hydrogen and helium. For the dust and gas 
opacity we set n{p, T) = K g p a T b , where Ko, a, b are 
constants (|Bell fc Linlll994D that depend on the species 
and the physical processes contributing to the opacity 
at each p and T, for example, ice mantle melting, the 
sublimation of dust, molecular and H~ contributions. 

Up to densities p ~ 10 _9 gcm -3 , the self-gravitating 
gas dynamics, energy equation and associated radiation 
transport are treated explicitly. Wherever a gravitation- 
ally bound condensation with p > 10 -9 gcm~ 3 is formed, 
it is presumed to be destined to form a protostar. Then, 
in orde r to avoid very sm all timesteps, it is replaced with 
a sink (B ate et al.lll995l ). i.e. a spherical region of radius 
1 AU. The sink interacts with the rest of the computa- 
tional domain only through its gravity and luminosity. 
Matter that subsequently flows into a sink, and is bound 
to it, is assimilated by the sink, and ultimately ends up 
in a protostar at the centre of the sink. 

The radiation from protostars that form during the 
simulation is taken into account by invoking a back- 
ground radiation field with temperature T BGR (r) (see 
Eq. [T]) that is a function of position relative to all the 
protostars present in the simulation (IsFa matcllos et all 
I2007bt [Stamatellos fc Whitworth! l2009bl ). 



r Bon M = (10 K) 4 + j jj— 



L n 



(2) 



where L n and r„, are the luminosity and position of the 
nth protostar. 
At all times, the luminosity of a protostar is given by 



r -i Mil 

1 M. 



JGM n M n 



(3) 



dui 
~dt 



4<7 SB (Tj 



4 

BGR 



rp4 



(1) 



where M n is the mass of the protostar, R n its radius, and 
M n is the accretion rate onto it. The first term on the 
righthand side is the intrinsic luminosity of the protostar 
(due to contraction and nuclear reactions in its interior) . 
The second term is the accretion luminosity. / = 0.75 is 
the fraction of the accretion energy that is radiated away 
at the photosphere of the protostar , rather than bein g 
expended driving jets and/or winds (Offn er et aLl l20091 
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We assum e R n = ZKp, is the ty pical radius of a young 
■ (|Palla fc Stahlerlll993t) . At the initial stages of 



protostar i[rana 6i STanienuaaji) . at, me miriai stages 01 
star formation the accretion luminosity dominates over 
the intrinsic luminosity of the protostar. 

3. A PHENOMENOLOGICAL TIME-DEPENDENT MODEL 
OF EPISODIC ACCRETION 

Simulations of self-gravitating hydrodynamics on the 
scale of molecular cloud cores (i.e. sizes from 10 4 to 
10 5 AU and masses from 0.3 to 10 M©) can achieve suffi- 
cient resolution to capture the formation of discs around 
young protostars, and the effects of GIs in the outer re- 
gions of such discs. However, such simulations cannot 
capture what happens in the inner disc region, where 
sinks are invoked. It is normally assumed that any mat- 
ter inside a sink flows instantaneously onto the central 
protostar. Here, we assume that the mass accreted into a 
sink is deposited on an inner accretion disc (IAD) inside 
the sink, where it piles up until it becomes hot enough 
that thermal ionisation couples the matter to the mag- 
netic field. At this juncture the MRI is activated, trans- 
porting angular momentum outwards, and thereby allow- 
ing the matter accumulated in the IAD to spiral inwards 
and onto the central protostar. Hence the sink mass is 
divided between the central protostar (M*) and its IAD 
(M 1AD ): 



(4) 



Matter is assumed to flow from the IAD onto the proto- 
star at a rate 



M* = 1O- 7 M yr- 1 



(5) 



where the first term on the righthand side is a low reg- 
ular accetion rate that obtains at all times, and the sec- 
ond term is a much higher accretion rate that obtains 
only when MRI acts to transport angular momentum in 
the IAD. The IAD is presumed to be a continuation of 
the much more extended accretion disc outside the sink, 

which is simulated explicitly. 

In detailed disc models (jzhu et alJl2009al fibL l2010bl |ah 
the matter in the IAD couples to the magnetic 
field, due to thermal ionisation, once the temperature 
reaches a thresho ld of T, JPT ~ 1400 K. U sing an a- 
parameterization (jShakura &; Sunvaevlll973D for the ef- 
fective viscosity delivered by the MRI. IZhu et al.l (|2010bl ) 
then estimate that the accretion rate during an outburst 
is 



(6) 



where a MRI (where ct„ PI is the effective Shakura- 
Sunyayev parameter dShakura &: Sunvaevll 19731) for MRI 
viscosity). iZhu et al.l pOlObt ) also find that the duration 
of an outburst is 



At MK1 ~ 0.25 kyr v (( } 



M* 
O.2M 

\ 1/9 



2/3 



10" 5 Moyr- 1 



(7) 



Here M IAD is the rate at which matter flows into the sink 
and onto the IAD. We assume that the critical tempera- 
ture for the MRI to be activated is reached when enough 



mass for an MRI-enabled outburst has been accumulated 
in the IAD, i.e. 



M lAO > M MR1 ~ M MR1 At 

MRI 



(8) 



Using Eqs. [Hand [7] we obtain 




We do not model in detail the thermal and ionisation 
balance in the IAD, since this is not the main concern of 
this paper. Instea d, based on observations and models of 
FU Or i-type stars (jHartmann fc Kenv on 1996; IZhu et al.l 
l2010bf) , we presume that rapid accretion onto the proto- 
star, at a rate given by 



M 

MRI 

MRI — 



exp 



(*-*o) 



,i o <i<t o +At MRI (10) 



is initiated as soon as M IAD exceeds M MRJ ; t Q is the time 
at which this occurs, and the outburst terminates at t Q + 

The time taken to accumulate (or re-accumulate) M MRI 
i.e. 



is At ACC - M MKI /M, 



: 13kyr 



0.2M f 



2/3 



M, 



-8/9 



10- 5 M Q yr- 1 



(11) 



This is much longer than the duration of the outburst 
(see Eq.0). 

With the above formulation the only free parameter 
is a MRI , which controls the strength and duration of 
the outburst; increasing a MRI makes for a more intense, 
shorter outburst. The mass delivered in each outburst, 
and the duration of the interval between outbursts, are 
essentially independent of ce MR1 , which is fortunate, in 
the sense that a MRI is rather uncertain. O bservations 
and s imulations suggest a MRI = 0.01 to 0.4 (Kin g et all 
12007th We have adopted a MRI =0.1 which is consis tent 
with observations of FU Ori events (Zhu et al. 2007). 

4. THE IMPORTANCE OF EPISODIC ACCRETION IN 
LOW-MASS STAR FORMATION 

To evaluate the consequences of episodic accretion for 
disc fragmentation and low-mass star formation , we per- 
form radiation hydrodynamic simulations of a collaps- 
ing turbulent molecular core, with properties chosen to 
match the observ ed properties of prestellar cores (e.g. 
And re et alj |2~000). The initial density profile is 



p{r) 



KERNEL 



1 + (r/R Y 

x ' \ / KERNEL / 



(12) 



Here p KERNEL = 3 x 10~ 18 gcm- 3 is the central den- 
sity, and -R KBRNBL = 5,000AU is the radius of the cen- 
tral region within which the density is approximately 
uniform. The outer envelope of the core extends to 
R CORE — 50, 000 AU, so the total core mass is M CORE = 
5.4 M . The gas is initially isothermal at T = 10 K, 
and hence the initial ratio of thermal to gravitational 
energy is a THERM = 0.3 . We impose an initial random, 
divergence-free, turbulent velocity field, with power spec- 
trum Pkdk oc fc~ 4 (ifc, to match the scaling laws observed 
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Fig. 1. — Evolution of the accretion disc around the primary protostar forming in a collapsing turbulent molecular cloud core, without 
radiative feedback from the protostar. The disc around the primary protostar increases in mass, becomes gravitationally unstable, and 
fragments to form 3 low-mass stars, 2 brown dwarfs and 2 planetary-mass objects. The colour encodes the logarithm of column density, in 
g cm -2 . 




Fig. 2. — Evolution of the accretion disc around the primary protostar with continuous accretion, continuous radiative feedback from the 
protostar. The disc grows in mass, but radiative feedback makes it so hot that it does not fragment. 
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Fig. 3. — Evolution of the accretion disc around the primary protostar with episodic accretion and episodic radiative feedback from the 
protostar. The disc becomes gravitationally unstable (first and second column) but fragmentation is damped by heating due to an accretion 
burst (third column); however, after this burst the disc cools sufficiently to undergo gravitational fragmentation; 2 low-mass stars form in 
the disc. 



in molecular clouds (|Larsonll"l981[ ). and amplitude such 
that a TURB = C/ TURB /|t/ GRAV | = 0.3. 

The core is assumed to evolve in isolation, i.e. any 
interactions with its environment are ignor ed. This ap- 
ears to be the case in Ophiuchus, where lAndre et ahl 
2007f ) conclude that individual cores do not have time 
to interact with one another before evolving into proto- 
stars. In denser star formation regions interactions be- 
tween cores may occur, but only very close encounters 
will influence disc fragmentation in scales less than a few 
hundred AU. Such encounters are ex pected to enhance 
disc fragmentation (|Thies et al.ll2010h . 

The molecular cloud core is represented by 10 6 SPH 
particles, so each SPH particle has mass 77i SPH ~ 5 x 
10 _6 Mq. The minimum resolvable mass is therefore 
M MIN ~ A/" NEIB ra SPH ~ 3 x 1O~ 4 M . This is much 
smaller than the thoretical minimum mass for star forma- 
tion (~ 3 x 1O~ 3 M0; usually refcrcd to as The Opacity 
Limit), and therefore all self-gravitating condensations 



should be reasonably well resolved. 

We perform three simulations, all with the same initial 
conditions, and differing only in their treatments of the 
luminosities of protostars. In all three simulations the 
initial collapse leads to the formation of a primary pro- 
tostar (i.e. a first sink) at t ~ 77 kyr, and this quickly ac- 
quires an extended accretion disc. The simulations only 
diverge after this juncture. 

4.1. No radiative feedback 

In the first simulation there is no radiative feedback 
from the protostar, and the gas in the disc is only heated 
by compression and by viscous dissipation in shocks. 
It is therefore cool enough to experience strong GI 
(Stamatellos fc Whitworthll2009a[ ). The resulting gravi- 
tational torques transport angular momentum outwards, 
allowing matter to spiral inwards and onto the primary 
protostar at a rate 10~ 5 to 10 -4 M Q yr _1 , but in the 
outer parts of the disc, at radii R > 50 AU, the GI also 
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Fig. 4. — The evolution of the primary protostar and its IAD. Red lines are used for the case with accretion but no radiative feedback, 
blue lines for the case with continuous accretion and feedback, and black lines for the case with episodic accretion and feedback. (A) shows 
the accretion rate into the sink and the accretion rate onto the primary protostar. For the first two cases these two quantities are the same 
(by definition), but for episodic accretion they are different. (B) shows the accretion luminosities for the cases with continuous (blue) and 
episodic (black) accretion; the luminosity is effectively zero for the case with no radiative feedback. (C) shows the mass of the primary 
protostar, which increases very steadily for the runs with no radiative feedback (red), and with continuous radiative feedback (blue). For 
the case with episodic radiative feedback (black), we plot the mass of the protostar and the mass of the IAD; the mass of the protostar 
only increases significantly during the brief accretion bursts. 



becomes highly non-linear locally, resulting in the forma- 
tion of seven secondary protostars, with masses ranging 
from 0.008 M Q to 0.24 M . This is illustrated in the se- 
quence of frames on Fig. 1. The accretion rate onto the 
primary protostar and the growth of its mass are shown 
on Fig. 4 (red lines). 



4.2. Continuous radiative feedback 

In the second simulation, we assume that the matter 
entering a sink is immediately accreted onto the proto- 
star at its centre. This results in an accretion luminosity 
which is typically 10 to 100 L©. It therefore heats the 
surrounding disc, and as a result GI saturates, gener- 
ating low-amplitude spiral waves that transport angular 
momentum outwards and allow matter to spiral into the 
protostar, but entirely suppressing disc fragmentation, 
at all radii. This is illustrated in the sequence of frames 
on Fig. 2. At all times the spiral structure in the disc is 
diffuse and has low amplitude, there is no fragmentation 
and therefore no secondary protostars form. The accre- 
tion rate onto the primary protostar, its resulting accre- 
tion luminosity, and the growth of its mass are shown on 
Fig. 4 (blue lines). 



4.3. Episodic radiative feedback 

In the third simulation, we use the semi-analytic, time- 
dependent model for angular momentum transport inside 
the sink that is described in Sect. [3] The luminosity of 
the protostar during the MRI-driven accretion bursts is 
very high, and stabilises the disc against fragmentation. 
However, as the mass of the protostar (M t ) builds up, the 
mass in the IAD that is required to activate MM (M MRI ) 
also goes up, and at the same time the accretion rate into 
the sink (M IAD ) goes down, so that the intervals between 
bursts become longer. During these intervals the lumi- 
nosity of the protostar is relatively low, hence the disc 
cools down. Within a few kyr of the formation of the pri- 
mary protostar, the interval between successive accretion 
bursts has increased to ~ 5 kyr, and this is sufficient time 
to allow the outer disc at R ~ 50 to 150 AU, to undergo 
gravitational fragmentation. This is illustrated on Fig. 3. 
In the first two frames, at 82 and 83 kyr, the disc is quite 
unstable and tries to fragment, but then at 84 kyr it is 
heated by an accretion outburst and stabilised. Once 
the outburst is over, the disc cools back down, and frag- 
ments to produce two low-mass hydrogen-burning secon- 
daries. (Other, statistically similar simulations also pro- 
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d uce brown dwarfs and even plane tary-mass objects, as 
in lStamatellos &: Whitworthll2009al) . The accretion rates 
onto the primary protostar and onto its IAD, its resulting 
accretion luminosity, and the growth of the mass of the 
protostar and its IAD are shown on Fig. 4 (black lines). 

5. DISCUSSION 

These results demonstrate the importance of episodic 
accretion in facilitating the formation of low-mass sec- 
ondary stars and brown dwarfs by disc fragmentation. 
The critical factor that determines whether disc fragmen- 
tation occurs is the time interval between successive ac- 
cretion bursts. In the phenomenological model presented 
here, w hich is based on the detailed model of iZhu et al.l 
(2010b]), this interval is given by Eq. [Til an d depends 
on the rate of inflow into the IAD, M IAD , and the mass 
of the primary protostar, M*. It is initially quite short, 
because M* is small and M IAD is large. However, over 
time M* increases, and in general M IAD decreases. By 
the time M* ~ 0.2 M Q and M IAD ~ 10 -5 M yr~\ the 
interval is ~ 13 kyr. This is sufficient time for the outer 
disc to fragment, since the condensation timescale for an 
unstable fragment at 50 to 150 AU is < 5 kyr. It follows 
that the relative abundance of very low-mass stars may 
be different in different star forming environments. In 
dense molecular cloud cores where the rate of accretion 
onto the IAD is expected to be high, due to the large 



reservoirs of gas available, the interval between accretion 
bursts may tend to be shorter, and hence the opportuni- 
ties for disc fragmentation less frequent. This should be 
reflected in a paucity of very low-mass stars and brown 
dwarfs in such regions, relative to more diffuse star form- 
ing regions. 

We conclude that if episodic accretion is a com- 
mon phenomenon among young protosta rs as observa- 
tional and theoretical evidence suggests (iHerbid 11977 



Dopita 1978: Reipurth 1989; Hartm ann fc Kenvonl 



1996 



iGreene et al.ll2008t iPeneva et al.ll2010D . it may limit the 
effect of the luminosity of a protostar on its environ- 
ment and create favourable conditions for disc fragmenta- 
tion to occur. As disc fragmentation is predominantly a 
mechanism for forming low-mass stars and brown-dwarfs, 
episodic accretion may then have a significant influence 
on the lower end of the initial mass function. 
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